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Abstract 

The kinetic theory of dilute gases to first order in the gradients yields linear relations between 
forces and fluxes. The heat flux for the relativistic gas has been shown to be related not only to the 
temperature gradient but also to the density gradient in the representation where number density, 
temperature and hydrodynamic velocity are the independent state variables. In this work we show 
the calculation of the corresponding transport coefficients from the full Boltzmann equation and 
compare the magnitude of the relativistic correction. 
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I. INTRODUCTION 



It has been recently shown that, when constitutive equations to first order in the gra- 
dients are introduced in the transport equations for a relativistic gas, the system presents 

n 

no intsability nor causality issues [lj. Moreover, the pathology identified by Hiscock and 
Lindblom |2j, which in part lead to ruling out first order theories, is due to the coupling of 
heat with acceleration. This coupling has a phenomenological origin and is in contradiction 
with the results obtained from relativistic kinetic theory. Due to this fact, the first order 
theories are being currently reexamined and proposed as solid frameworks from where one 
could extract the physics of high temperature systems present both in astrophysical and 
experimental scenarios 

In this paper, the heat flux constitutive equation is obtained together with the associated 
transport coefficients using the full collision kernel in Boltzmann equation to first order in 
the gradients using a representation where the density, hydrodynamic velocity and temper- 
ature are the independent state variables. The heat flux in this scenario is coupled with 
temperature and density gradients and two transport coefficients are identified. The explicit 
form of such coeficients are obtained for a constant scattering cross section model and the 
results are shown to be consistent with the ones obtained by other authors using a different 
representation only in the comoving frame. However, for an arbitrary observer, the stress 
energy tensor includes Lorentz transformation factors The form of such tensor, from 
which the transport equations are to be extracted, is briefly discussed. 

The present work is divided as follows. In Section II the theoretical framework that 
sustains the calculation is presented. In Section III the mathematical problem is set up 
as two separate integro-differential equations whose solutions are formulated as expansions 
in orthogonal polynomials from which the general form for the transport coefficients is 
obtained. A constant scattering section model is assumed in Section IV in order to calculate 
the collision integrals and compare the transport coefficients as functions of z. The discussion 
of the results, including their comparison with the solution obtained in Ref. 4] as well as 
final remarks regarding the form of the stress-energy tensor are included in Section V. 
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II. SPECIAL RELATIVISTIC BOLTZMANN EQUATION 



Consider a neutral, dilute, single component, non-degenerate relativistic fluid in the ab- 
sence of external fields. Additionally, assume that such a gas is characterized by values of 
the relativistic parameter z, defined as the ratio of the thermal energy to the rest energy a 
single particle, close to one. Much higher values of this parameter correspond to very high 
temperature, ultra-relativistic gases for which the neutrality of the particles here assumed 
may not be appropriate. The non-relativistic gas corresponds to the limit where z tends to 
zero. Thus, for the system here considered, the space-time is given by a Minkowski metric 
with a + + H — signature for which the position and velocity four-vectors are given by 

x v = (x, ct) v v = 7( tt ) (w, c) (1) 

where c is the speed of light and 7( to ) = (1 — w 2 /c 2 ) 

The distribution function / is an invariant such that / (x u , v v ) d 3 xd 3 v , also an invariant, 
is the number of particles contained in a volume in phase space. The relativistic Boltzmann 
equation for the evolution of / reads |4-6|, 

V a f,a = J (f, f) (2) 

where the contraction on the left side corresponds to a total proper-time derivative. That 
is, in the molecule's rest frame where w = 

v a fa = ±d_L = dJ_ (3) 
c dr dr 

The right hand side of Eq. (j2j) is given by 

J(ff)= [ [ {f'fx'-ffx}Fa{n)dmv{ (4) 

n 

where dv^ = d 3 v\/vf and T is the invariant flux given by |4| 

T = ^2«M = \\l (^ia) 2 ~ ci = \^ (7(«07(iui) (w-wi- c 2 )) 2 - c 4 (5) 

which reduces to the relative velocity in the non-relativistic limit. The solution to the ho- 
mogeneous Boltzmann equation is determined by J (//') = together with the requirement 
of consistency with the local equilibrium assumption. That is, thermodynamic equilibrium 
is locally assumed and thus the state variables are given by 

N v = / f w v u dv* (6) 



which are the particle flux and equilibrium stress-energy tensor respectively. The thermody- 
namical, local equilibrium, variables for the system can then be extracted from the previous 
tensors as 

n = — ( 8 ) 

ne = M^T^. (9) 
c 2 

It has been recently shown that the heat flux can be defined, as in the non-relativistic case, 
as the average of the peculiar or chaotic kinetic energy Thus, since this velocity is the 
one measured by an observer locally comoving with the fluid element, the calculation will 
be performed in a comoving frame. Therefore, the fluid's hydrodynamic four velocity has 
only a temporal component: 

U u = (o, c) (10) 

The dissipative terms in the stress-energy tensor, as defined in Ref. Q from a standard 
tensor decomposition include a four- vector that in an arbitrary frame can be calculated 

as 

^ = c 2 L^q u (11) 

where q v is the heat flux, calculated in the comoving frame, and a Lorentz transfromation 
between the laboratory and each local equilibrium element of the fluid. This idea was firstly 
set forward, for the equilibrium quantities, by S. Weinberg [8]. This result was obtained 
by introducing such transformation to relate chaotic and molecular velocities and shows 
that the heat flux can be consistently defined only in the comoving frame. Because of that, 
from now on we will consider the hydrodynamic four-velocity as given by Eq. ( TTOl) and the 
molecular velocity v M will correspond to the chaotic velocity. 

In order to solve Eq. (|2j) the standard Chapman-Enskog method will be used J4, {J. Thus, 
the solution is approximated by 

/ = / (0) (l+^M) (12) 

where is the Jiittner equilibrium distribution function which, in the comoving frame, 
reads |l0| 









' rwyv 4 ' \ 




t mv u J 



Here z = kT/mc 2 is the relativistic parameter, where T is the local temperature, k the 
Boltzmann constant, and K n is the n-th modified Bessel function of the second kind. The 
solubility conditions imposed on <f>(v) are given by 

dv* = (14) 

which amounts to restrict the local state variables to be defined through the local equilibrium 
state. The proposed solution given in Eq. (TI~2|) is substituted in Eq. ([2]). Considering the 
deviation from the local equilibrium state <fi(v^) to be a first order quantity, one obtains a 
linearized first order Boltzmann equation which can be written as 

^ / (O) = / (O) C(0K)) (15) 

where the linearized collision kernel is given by 

C {</>) = [ [ {0' x + 0' - 0x - 0} /} 0) To (O) dQdv* v (16) 



The general solution to equation (fT5|) is given by the sum of the homogeneous solution plus 
a particular solution, = 0# + (pp. The homogeneous solution is obtained as a linear 
combination of the collision invariants 

1=0. (17) 
m^v J 

Existence of the particular solution is guaranteed by imposing an orthogonality condition 
on the homogeneous solution and the inhomogeneous equation namely, 



rwyv 4. 



v a f^dv* = 0, (18) 



quations ( JT8|) are the relativistic Euler equations obtained through the equilibrium solution 
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12j. In the absence of external forces, the left hand side of the relativistic Boltzmann 



equation is written 

af a fdf w ,df<f» A f . 

«/,«■=« + ~df~ > a + ft^fo ) ■ ( 19 ) 

The next step consists in substituting the derivatives of the Jiittner function and using the 
Euler equations to write the time derivatives in terms of the gradients. Such equations 



constitute a closed set for the state variables. At this point an appropriate representation 
needs to be chosen and thus we consider n, T and U v as the set of state variables. The Euler 
equations are written as 

h = -nu^ (20) 



K 3 V n T 



where () = u u (). u . Here the gradient of the hydrostatic pressure has been written in terms 
of the gradients of the number density and temperature by using an ideal gas equation of 
state which can be easily shown to hold for dilute special relativistic gases. We have also 
used the relation 

p _MY) (23) 



^ (f + *) K 3 (I) 

which can be verified by calculating e and p from the local equilibrium distribution function. 
After a somewhat tedieous but straightforward algebraic manupulation one can write Eq. 
f fT5|) as follows 

^?l^fl-7^V??fl + 2 -7§4fl-^IHCW (24) 



n 



MS) T \ z MS) *MS) 



where the term proportional to the hydrodynamic velocity gradient does not arise since the 
calculations are performed in a comoving frame. 
The solution to Eq. (I24p is given by 

1 M n 



<f> = A{i) v?h%^ + B (7) + a + a u v v . (25) 



The first two terms are the particular solution and the last two terms correspond to the 
solution of the homogeneous equation. The solubility conditions are thus written as 

A^yh^ + B^yhf^- + a + a v v v \ ifif^dv* = (26) 

where if) = mv tJ, ,m'y 2 . These conditions imply, as shown in Appendix A, that the constant 
a vanishes and cup is proportional to both h^n ya and h^T^ a so that Eq. (1251) reads 



4 = Abyufe + Bfryh^ (27) 

1 H n 
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In Eq. ([27)1 . since n Q and T Q are considered independent forces, A (7) and £> (7) are subject 
to the constraints 

j A^Wf^dv* = 0, (28) 

J B(^ 2 u 2 f^dv* = 0. (29) 

To take full advantage of the fact that the unknowns A and B are functions of 7, we perform 
all integrals in such variable using the relation 

dv* = Anc 3 ^-f 2 - ldj (30) 

which is obtained in Appendix B. 



III. EXPANSION IN ORTHOGONAL POLYNOMIALS 



^ 1 - ^IFlfi = C(B(7)^) (32) 



By substituting Eq. (]2Tj) in Eq. (|24p . the mathematical problem set up in the previous 
section yields two independent integral equations given by 

^{ l+ " _7 tf -M)} = ^ yhf) (31) 

and 

Mi) 

(;) 

subject to the contraints given by Eqs. (128]) and (129]) respectively. The unkown coefficients 

A and i3 are written in terms of orthogonal polynomials in 7 

00 

-4(7) = ^a n £ n ( 7 ) (33) 

n=0 

00 

8(7) = E ^ 34 ) 

n=0 

which satisfy the orthogonality condition 

£„(7)£ m (7)p(7)d7 = 5 nm , (35) 

where the weight function ^(7) = exp~?(7 2 — l) 3 / 2 . In the case where the hydrodynamic 
velocity is the one given in Eq. (fTUl) . these polyomials, the first two of which are obtained 
in Appendix C, are related to Kelly's set (4, -R3 by the relation 



R\ = \ l3K 2 ( - U£ n ( 7 ) (36) 



2 \/ \ Z 



In terms of the polynomials C, the subsidiary conditions can be written as 

oo „ 

Y^Cln / C n (j) P (^ = 0, (37) 
n=0 ^ 

oo „ 

J^K / £ n ( 7 )p(7)^7 = 0, (38) 



n=0 

Since £0(7) is constant, we have that a = b Q = and thus 



^( 7 ) = ^a n £ n ( 7 ), (39) 



n=l 

OO 



B(j) = J2 b " C M- (40) 



n=l 



The heat flux in the Chapman-Enskog approximation, as clearly stated in Ref. [3|, is given 
by the average of the chaotic kinetic energy flux, a definition that encompasses the physical 
conception of heat since the early developments of kinetic theory Since in this work 

the molecular and chaotic velocities coincide, we can write 

q » = mc 2 h^ [ jv u f {1) d*v (41) 



or, subtituting Eq. 



where 



nm 

q ~ AnczKoi- 



J (a) + 7 (6) 



(42) 



Jfa) = Ktfrf f Jv^A^e-^v (43) 

% = Ktfr* J 7 ^VS( T )e-^V (44) 

Notice that in both integrals only the u, /3 = 1, 2, 3 terms survive and from them, all v 7^ (5 
ones also vanish because the integrands are odd in the three-velocity. Thus, introducing 
Eqs. ( 1391) and (I40|) . the integrals read 

.5 



1(a) = —^^f E fl « / 7A»(7)p(7)d7 (45) 

n=l 

% = V E 6 " / T>Cn(7)p(7)rf7- (46) 

n=l 

As shown in Appendix C, we can write 7 = Co£o( 7 ) + c i^i( 7 )) with ci = ^/3g (z)z where 

9W = 5,A- 3 (i) + A- 2 (i)-M| (47) 



Using Eqs. (I45H47P in the heat flux given by Eq. (j4"2"|) we obtain that 



1 n 



(48) 



where the coefficients appearing in Eq. (|4~8]) are defined as 



nmc A ci , . 

= -^KJT) bl - (50) 



IV. SOLUTION OF THE INTEGRAL EQUATIONS 

The coefficients c^and b\, in terms of which the coefficients in Eqs. (149 D and (}5T)j) are 
given, have to be obtained from the solution of the integral equations (l3~Tj) and (132]) . In this 
section we outline such calculation to a first approximation. The variational method used 
is the standard one as described in detail in Ref. |16] . 

First we notice that the integral equations (I31H32P can be written as 



" VM?)? ~ i) f^^^ = ^E«nC(£„( 7 )^/ (0) ^ 



(51) 



- -^^^1(7) = ^^6nC(£ n ( 7 )^/ (0) ) (52) 

such that both have a similar structure. Indeed, since the dependence on 7 on both is 
the same, the procedure only needs to be carried out once for one of the equations and 
the solution for the other one can be readily inferred by adjusting the dependence on the 
parameter z. This similarity is consistent with the calculation in Ref. [4] where only one 
integral equation needs to be solved for the coefficient of a generalized thermodynamic force 
which includes contributions from VT and Vp in a single term. 

Following the presciption mentioned above, we will only deal with Eq. (152]) . Multiplying 
it by LrrSc^Vyh 1 ^ and integrating on both sides 



- £^ f h^vyC.C^v = h^b n J C m v v C{Cy f^)d% (53) 

where we have omitted the z and 7 dependences to short notation. For the integral on left 
hand side, using that e~~ *h v pV v vP = 47rc'p('y)d'y, we have 

j e-ltfpvyCxWCnWcrv = 47rc 5 5 lm (54) 



and thus, denning the collision brakett in the standard way 



[G,H] 



1 



G a 



H al + H a - H al - H a f^f^Ta (0) dQdvldv 



(55) 



in this equation G a = G a (vp) and we have used H al to denote H a {ya-^. The integral 
equations can be written as 



n=l 



^3 



z / n 



Mi) 



hp^bn [£ m {l)v v , C n {i)v p ] = c 2 ^f g [f[ 5i m 



(56) 



(57) 



n=l 



' ^3 (J) 

By following the standard variational method [16], the first approximation for a\ and b\ can 
be shown to be given by 



6l = c2 ^(| { ^ [^(7)^>A(7)^]} _1 



(58) 



(59) 



Thus, in order to calculate the coefficients a\ and b± to this level of approximation only one 
collision integral needs to be calculated namely, [£1(7)1?,,, £1(7)1^] . In order to calculate 
such brakett, the well known identity 



[G a , Hp] 



1 

An 2 



G a h/ + G a < - G a K) - G a «) 



(60) 



Hp ( < ) + ^ ( ) - ^ «) - Hp (v") f^f^Fa (A) dfid^W 



will be used. Using also the momentum conservation law for collisions and after several 
algebraic steps one can show that 

Xi-X 2 



h v p [A(7K,A(7)<1 
where the integrals Ii and I2 are defined as 



3g (z) z 2 n 2 m 4: c b 



(61) 



-7- 47 
±1 = m c 



J 2 = _ m 4 c 5 




7 



hX + (7 2 )' - (A - 7 2 1 fMfPTatmvlfo 




V [( 1V p)[ + (jvp)' - (jvp), - 1V p\ fWfl 0) Fadndvldv* 



(62) 



(63) 
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Thus, the general expressions for the coefficients Lt and L n are given by 

^_ 3mw^(S) A ( g(z) V (64) 



3£; 2 T 2 n 2 m 3 c 7 (#(z)) 2 



(65) 



Xx-2b ^Q)^(i)- 
In order to evaluate the integrals in T\ and X 2 , a particular collision model needs to be 
proposed. In the next section the simplest collision model, namely a constant cross section, 
will be assumed in order to obtain expression for the coefficients and asess their relative 
magnitude. 

V. CONSTANT SCATTERING CROSS SECTION 

The simplest model that one might consider in order to calculate collision integrals con- 
sists in assuming a constant cross section. The details of the calculations for such model can 
be found in Ref. |4J]. Here, we will only quote the final results for the integrals in Eqs. (162]) 
and fl63D 

647m 2 fc 6 T 6 CT / (2\ 1 /2\\ , . 

X i = : 7777 77772 2 ^ - + -K 3 - (66) 



{Mm 



647rn 2 fc 6 T 6 cT /4 /2\ 1 /2\\ , . 

" reV (g 2 (^)) 2 V* W + ? W/ <67> 

where a is the constant scattering cross section. Subtituting these expressions in Eqs. ([6 
and (1651) one obtains 



3c*r/ jf a (i) \ (Mi) -10(91) 



(68) 



(69) 



64- *K 3 ® + + 2) K 2 ®) 

which are to some extent approximations to the coefficients in appearing in the heat flux con- 
stitute equation in the n, T, U u representation. In the limit of small relativistic paramenter 
z one recovers the non-relativistic values 



J TNR 



I- Y Q Z + (70) 
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L n ~ -L TNR ^z- — z 2 + ...j (71) 

where Ltnr = (75mc 3 /256y/na)z 3 ^ 2 is the usual non-relativistic thermal conductivity for 
hard spheres divided by T. Notice that indeed, as z goes to zero Lt —> Ltnr and L n — > 0. 
In Fig. [I] the magnitudes of both coefficients, normalized to Ltnr are shown as functions of 
z. 



1.2 r 




0.2 - ^ 

' 

/ 

J 

1 2 3 4 5 



Figure 1: The magnitudes of both coefficients, normalized to Ltnr a re shown as functions of z. 



VI. DISCUSSION AND FINAL REMARKS 



In the previous section the transport coefficients involved in the constitutive equation for 
the heat flux in terms of the forces VT and Vn have been obtained. The calculation has 
been performed in the fluid's local comoving frame, that is in terms of the chaotic or pecular 
velocity. The Chapman- Enskog solution method of Boltzmann's relativistic equation in this 
representation leads to two independent integral equations. Once the solution is written in 
terms of a suitable orthogonal set of polynomials, both equations have the same structure 
and the solution is finally obtained using the standard variational method. By assuming a 
very simple constant cross section model to calculate the collision integrals, the coefficients 
can be studied. The results are shown in Fig. 1. It can be clearly seen that for small values 
of z the non-relativistic limit is verified. For larger values of the relativistic parameter 
the coefficient L n increases and has the same order of magnitude of Lt- That is, if the 
gradients of the system are of the same order of magnitude the relativistic effect q oc Vn is 
as important as the usual Fourier effect. 
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In previous work, Cercigniani & Kremer obtained the constitutive equation for the heat 
flux in a relativistic neutral gas in a different way. The definition they used for the heat flux 
is the one obtained from the phenomenology using Eckart's decomposition in an arbitary 
frame. However, in the comoving frame both definitions, the one given here in Eq. (|4l|) and 
the one they consider namely, 



(72) 



where p a is the four-momentum, are the same. However, after the calculation of the heat 
flux from the Chapman-Enskog solution of the Boltzman equation, in Ref. 4| the heat flux 
is obtained in terms of a relativistic thermal force which includes both the Fourier term as 
well as the relativistic correction which they leave in terms of the gradient of the hydrostatic 
pressure. Their result is 



T 
nh f 



where 



and A is given by 



e+ t = P I<3 (j) 
n nz K 2 Q) 



» 5 (l + 5f?(i)-l0(i) : 



h - c 2 h 



(73) 



(74) 



(75) 



with the integrals Ii and I 2 defined as 



I\ =UaUpUyUsX. 



I 2 =U a U^x 

P> P {(P»; + (PW - (P» x -P» f^f^Fadn^^- 



pl p^ 

The invariant flux in the expressions above is related with the one in Eq. <^ by F — m 2 cJ-" 
and also the distribution function in Ref. [4] differs by a factor from the one used in this 
work: f r = 4rf. 

In order to compare Eq. (173]) with our results the hydrostatic pressure can be written in 
terms of Vn and VT by means of the ideal gas equation of state. Also the hydrodynamic 
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velocity has to be substituted by the one given in Eq. (ITU1) . which accounts to writting 
Cercigniani's result [4] in the comoving frame. After these changes have been introduced, 
the integration variable changed to 7 and the different signature is taken into account, one 
can readily verify that I\ = T\ and I2 = c~ 2 X 2 and that the constitutive equation in Eq. 
( I73p can be written as 



h a0 



Ml)) T K 3 (l)\n 



(76) 




which has the same structure of Eq. (|48)) . Subtituting the value for A obtained in Ref. j4j 
which can be written as 

A = _^__ ( ^_, (7?) 

yields 

A7Y1 2 Mi)\_(Mli) !\ 3nWmV / g(z) 

XT v " z md ) x^m. 1 (78) 

AT M$ = 3n 2 fc 2 T 2 mV (g {z)f 

Mi) J i- c2/2 MS) Mi) 

which are precisely the coefficients in Eqs. (l64l) and (155]) . Thus, one concludes that both 
calculations are consistent in the comoving frame. 

As mentioned above, the constitutive equation obtained in this work is equivalent to the 
one obtained in Ref. |4| only in the comoving frame since in that case the definitions for 

n 

the heat flux are identical. In this work we follow the ideas set forth in Ref. [3| where 
the heat flux is a local quantity which only makes sense in the comoving frame where the 
molecular velocity coincides with the peculiar velocity. However, the term appearing in the 
stress-energy tensor and that will ultimately impact the transport equations is in our case 
given by = c 2 L^q v with q v given by Eq. (1411) while in the traditional calculations what 
is obtained is r M = q v with q v given by Eq. (1721) . This difference is discussed in Ref. [3!] and 
its implications will be explored in future work. 
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Appendix A 

In this appendix we show that Eq. (127)1 is valid in view of the subsidiary conditions given 
in Eq. (l2ljj) . Equation (I26p can be written as 

■A^y Y + B (l)v ir ^ + a + &d + «4f 4 j ijf {0) dv* = 0, (80) 

where we have separated spatial and temporal terms in the contraction a v v v since they have 
opposite parity. 

For ip = v A = 7 c the first, second and forth terms yield odd integrands and thus 

J (a + & iV A ) jf^dv* = (81) 

For ijj — v with I — 1, 2, 3 the third and fifth terms do not contribute because of their odd 
parity and we obtain 

/ (a( T )^ + B(7)^ + = 0, (82) 

which vanishes for £ ^ k. Thus 

J U 7 ^ + B^ + vWfVdv* = ^ J + B^ + a e ) jWf^dv*, (83) 

and thus the condition in Eq. (1821 reduces to 

/ U(j)^ + S( 7 )^ + fi«) jWf^dv* = 0. (84) 
Finally, for ip = 7 2 the parity in the terms is the same as for ip = v 4 such that 

j (a + a 4 v 4 ) 7 2 / (0 W = (85) 

From equation (182]) we have that 

_ T, k / A^W^f^dv* n, k f B^tu'f^dv* 

ak T f^uj 2 f(°Uv* n J^u 2 f^dv* ' 1 J 

such that we can redefine 

Al)v k7 ^- + B(j)v k ^ + a k v k — ► a( 7 )^ + 6(7)^— (87) 
in in 



and now the subsidiary condition (182|) reads 

n 7 V/ (0) ^* = 0. (88) 



A( 7 )^ + £( 7 )^ 
1 n 
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Equations ( IHTj) and (185|) can be written as an homogeneous system for (a, d^): 

ag n + d^gvi = 0, 

ag 2 i + = 0, (89) 

where 

9n = J lf (0) dv*, g 12 = g 2l = j <fcf®dv*, g 22 = J -f<??®dv\ (90) 

such that the determinant does not vanish and the solution is the trival one a = d\ = 0. 
Putting toghether these two results, we can conclude that the proposed solution consistent 
with the subsiadiary conditions is the one given by Eq. fl27j) where the condition in Eq. (1821) 
still needs to be enforced. Considering both forces and n & as independent forces, this 
requirement is written as two separate conditions in Eqs. ( 1281) and ( |29l . 



Appendix B 



In this appendix, the identity dv* = Atcc 3 ^^ 2 — ld'j, which is used in several parts of the 



work, will be obtained. The invariant volume element in velocity space is given by 



HQ 



or, in terms of the three- velocity w 



dv* = c d 4 = d jl (91) 



d 3 w 

dv* = det [J] — (92) 

7 



where the Jacobian matrix has components = dvi/dwj which can be shown to be given 
by 

^=7fe+7 2 ^) (93) 



Using the identity det [<% + AiBj] = 1 + A i B i 



,w 2 



det[J]= 7 3 (l + 7 2 ^) = 7 5 (94) 



si 

obtain 



and thus, introducing spherical coordinates for w and that w 2 dw = ■% \/j 2 — 1g?7 we finally 



dv* = 47rcV7 2 - ld l (95) 
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Appendix C 



The orthogonal polynomials are obtained using the standard Gram-Schmidt procedure 



17j. The proposed polynomials are C n {^) = a 0n + a in 7 + a 2 „7 2 + • ■ ■ + a n nl n which are 



required to satisfy the orthonormality condition 



71 



(96) 



where p (7) = exp I ('y 2 — l) 3 / 2 . For n = we have Coil) = a oo an d the orthonormality 
condition ( 196]) yields 



A)(7) 



For n = 1 we have £1(7) = aoi + a n7 an d equation ( |96|) yields 



(97) 



a i 



Mi) 1 

1 



I 1 ) •'-••(' 



On 



v 7 ^ 



z J \ : 



Mil! 



Mi) 



(98) 
(99) 



if we define g (z) = 5zK 3 (±) + K 2 (I) - , ( _, 

A(7) ' 



^3 f 1 ) 

TY+7 



(100) 



\/ 3 # 

We only need these two polynomials in order to calculate the heat flux. Additionally, the 
coefficients in 7 = CqCq + C\L\ are introduced in Sect. III. A simple way of obtaining them 
is to solve Eq. f llOOp for 7 and use Eq. fl97|) which yields 



7 = y /3g(z)zC 1 ^) + y/3z , 3{ ~ J Co 



Mi) 



(101) 



from which c = \[3z ^ 3 ^\ and c\ = ^3g (z)z. 
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